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INTRODUCTION 



Metal foams are a new class of materials with low densities and novel physical, mechanical, thermal, 
electrical, and acoustic properties. They offer potential in fuel cell systems for thermal management, 
catalyst support, filtration, steam generation, mixing, energy absorption, acoustic absorption, and 
filters as well as applications in the cell stack assembly as bi-polar plates and gas diffusion. Metal 
foams also hold particular promise in applications where several of these features can be exploited 
simultaneously. 

Porvair Fuel Cell Technology produces a variet}^ of porous structures (Figures 1&2). The base 
materials used to make them can be most any powder metal (or ceramic), the void fraction and pore 
size can be controlled independentiy of one another, and virtually any shape can be manufactured. 
The pore sizes can range from about 0.2 inches to 10 microns and the void fraction can range from 
about 97% to 80% (3% relative density to 20% relative density). This wide range of base materials 
and properties guarantees a product that can fit in most any application. 



Figurel: FeCrAlY Metal Foams in 

Catalyst Substrate and Heat 
Exchange Applications 



Figure2: Copper, Brass, FeCfAly, 

Nickel, Titanium and Stainless 
Steel Metal Foam 




Controlling the material properties is very important for many different reasons. For instance, a very 
small pore size has a larger surface area for catalytic applications than a larjge pore size, but comes 
with the cost of a higher pressure drop. A foam structure with a high relative density has better 
conductive and convective heat transfer properties as well as greater strength than foam with a low 
relative density. But an increase in relative density means a heavier component with a slower time 
response in thermal applications. 

In order to use metal foams effectively, the interaction between the physical characteristics (relative 
density, pore size, material type) and the material properties (strength, pressure drop, surface area) 
must be understood. This paper attempts to present the properties in a way that facilitates 
comparison with other materials and structures, summarizes the guidelines for design with metal 
foams, and illustrates how they might be used in fuel cell applications. 



PRESSURE DROP IN IRON BASED METAL FOAMS 



Experimental Setup and Procedure 

The pressure drop tests apparatus is shown schematically in Figure 3, it consists of a flow meter, 
pressure gauges, and hand operated valves. Air flows through die inlet pipe which is four inches in 
diameter and 30 inches long. The inlet pipe terminates at die test section, where the specimen is 
loaded. The specimen is placed into a gasket, to ensure die air flows in the axial direction, and is 
secured widi a pneumatic clamping system. The two valves shown m Figure 3 are adjusted until the 
desired flow rate is obtained. The inlet and exit pressures are noted as weD as the flow rate. Based 
on die results from die measurements, die pressure drop as a function of velocit}^ is obtained. 

Figure 3: Test Apparatus Schematic 
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Equation Development 

Darcy may have been the first to examine flow through porous media. Darcy [1] described his 
empirical studies on steady-state filtration demonstrating that the water volume passing through a 
sand layer is proportional to die pressure-drop across the length of die layer, leading to the equation 



u' 
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where a is the hydraulic conductivit)^ of the porous medium, U* is the channel cross-section average 
speed, and AP' is the hydrostatic pressure difference across the layer of lengdi L". 

Several steps have been taken to generalize Eq. (1). Hazen [2] unveiled the fluid viscosity effect on 
die original flow equation proposed by Darcy, Eq. (1). By altering die temperature of the fluid 
before entering the medium, he observed the influence of the temperature on die hydraulic 
conductivity calculated using Eq. (1). He proposed a modification to die original equation of die 
form 
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where T is die water temperature, measured in degrees Fahrenheit, and aso is die reference hydraulic 
conductivity value measured when the water is 50°F. 



Inspection of Eq. (2) reveals that the correction 60/(T*+10) mimics the water viscosity dependence 
on temperature, fX(T*)/|Ll(500F), to less that one percent. Therefore, it is reasonable to rewrite Eq. (2) 
as 



ti(T ) L 

where K=[|U-(500F) ttso], known as the specific permeability, is a parameter supposedly independent 
of fluid properties. 

Dupuit [3] was the first scientist to use theoretical principles to interpret Darcy's equation [4]. A 
consequence of Dupuit's study is the inclusion of a quadratic velocity term added to Eq. (3) modeling 
the resistive form drag imposed on the fluid by the solid permeable medium 

^^mu-^cpu-'. (4) 

L K 

where the coefficient C is the form coefficient of the porous layer. Eq. (4) is referred to as the 
Hazen-Dupuit-Darcy equation [4]. 

In an internal technical memo [5], J.A.F. Hill reported the development of an equation to predict the 
pressure-drop across a porous-metal filter as a function of fluid flow. The pressure-drop data was 
gathered by flowing air though a porous metal sheet. The origin of the equation can be found in 
publications by Forscheimer and Reynolds, In the basic equation, 

m is the mass flow rate, and all other variables are as mentioned before with a slightiy different 
interpretation. The first constant, K, in Eq. 5 is called the inertial resistance coefficient and is 
analogous to the specific permeability. The second constant, C, is called the viscous resistance 
coefficient and is analogous to the form coefficient. At low flow rates the contribution of the C term 
is negligible, but as die flow rate increases, the contribution of C is significant. At very high flows C 
dominates the pressure drop relationship. 

By rearranging Eq. (5) and using the equation of state for ideal gas, P=pRT, we get 
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Integrating Eq. (6) across the bed length, L*, and assuming constant fluid viscosity gives 



^ = ^^^=^m + 2RTCm\ (7) 
K 



Rearranging Eq. (7) yields 
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but according to M.J. Rosso Jr. et aJ [6], Eq. (8) takes the form 

L K 



(9) 



after "suitable fitting to measured results". Rosso et al indicates that K is a funcdon of the system 
geometry as well as die characterisdcs of the porous material. Unlike Eq. 4, Eq. 9 implies that AP^ is 
proportional to CiU'-i-CaU'- rather than AP. 

Rosso et al also indicates diat K should be proportional to 1/Dp- and Prer/(l-prci)^. Where Dp is the 
particle diameter and prd is the relative density (1 minus the void fraction) of the porous material. It 
should be noted that the experiments by Rosso et al were conducted on porous metal sheets not metal 
foam structures. The particle diameter, Dp, does not play a significant role in the pressure drop of 
foam structures; the pressure drop is primarily driven by the open pore size (cell size). 

In the book, Metai Foams: A Design Guide Ashby et al [6] show the pressure drop as 

a a 
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where a is the metal foam cell size, Va is the fluid kinematic viscosity, Pa is the fluid density, a is the 
fluid thermal diffusivity, U is the free steam velocity, and d is the ligament diameter (or strut 
diameter). The exponent m and the coefficient ^ are determined experimentally. 

For the purposes of this study, Eq. 10 was rewritten with ppi number and relative densit)^ We know 
that the ligament diameter is inversely proportional to the ppi number and the relative densit)^ Also, 
the metal foam cell size, a, is inversely proportional to the ppi number and relative density. 
Therefore Eq. 10 is rewritten as 
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Equations 4, 9, and 1 1 are utilized and the results are compared in order to shown the dependence of 
pressure drop on the free steam fluid velocit>% the relative density, and the ppi number. 



Experimental Results For Constant Inlet and Exit Diameter 

Scrutiny of Eq. 4 

Nearly 180 runs were conducted with stainless steel samples of different diameter, thickness, inlet 
velocity, ppi, and reladve density. Equations 4, 9, and 11 are utilized and the results , are compared to 
show the best fit for the data. Figure 4 shows the interaction between the pressure drop and the inlet 
velocity for 5 ppi stainless steel foam with different relative densities. The same plots were created 
for 10, 20, 30, 40, 50, 60, and 80 ppi foams with similar results and are located in Appendbc I. The 




solid lines are a quadratic least squares fit of the data and show an average coefficient of 
determination value of 99.6%, indicating that the model relating pressure drop to the velocity can 
account for 99.6% of the variation present in the sample of pressure drop values. Moreover, the 
qualitative results from the curve fit shown that as the ppi number and relative density increase so do 
the coefficients 1/K and C in Eq. 4. The quantitative results of die curve fit in conjunction with xht 
qualitative results tend to validate the use of Eq. 4. 

The values of ^l/K and pC were computed at various densities for each ppi number. The values are 
tabulated in Appendix II and show an increase in the two constants as the relative density is 
increased and as the ppi number is increased. The constants were then fit to the relative densit}^ at 
each ppi number. Figure 5 shows the relationship between the relative density and |^/K, as well as 
the relative density and pC at 5 ppi. 
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Figure 5: mu/K and rho C as a Function of Relative Density 




A linear dependence of Ji/K and pC on the relative density was initially assumed. This lead to the 
equations: 

|: = «iP«/ (12) 

and 

CP=CC2Prei^ (13) 



where tti and az are functions of the ppi number. The coefficients tti and tta were then fit to the ppi 
number with the results shown in Figure 6 and Figure 7. 



Figure 6: CL\ as a Function of PPI Number Figure 7: CLz as a Function of PPI Number 
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The coefficient ai was found to have an exponential dependence on the ppi number while the 
coefficient a2 was found to have a linear dependence on the ppi number. Finally, an equation was 
developed to quantify the coefficients of Eq. 4. They are as follows: 

■^ = 19.5734p„/°'"<'""' (14) 
K 

Cp = 2.9359Sp,,,ippi) (15) 

Figure 8 shows the relationship between the actual pressure drop and the model generated from Eq. 
4, Eq. 14, and Eq. 15 (Model I). The chart indicates a coefficient of determination of 91.2% 
implying that the model accounts for approximately 91% of the variation in the actual data. A 
slightiy more complex adaptation of Eq. 14 and Eq. 15 was explored yielding a coefficient of 
determination of 96.2%. The adaptation is as follows: occurrences 

i^- = 13.6173p^re°°^"<'""""* (16) 
K 

Cp = 2M6p'J'ppi''' (17) 

Figure 9 shows the relationship between actual pressure drop and the model generated from Eq. 4, 
Eq. 16, and Eq. 17 (Model II). 
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Figure 8: Model I (Eq 4, 14, & IS) vs Actual 



Figure 9: Model II (Eq 4, 16, & 17) vs Actual 



250000 
? 200000 

r 150000 

o 

I 100000 



50000 





7 . 



















100000 200000 300000 
dP/L Actual (Pa/m) 



250000 
200000 

1 

r 150000 

o 

I 100000 



% 50000 




100000 200000 300000 
dP/L Actual (Paint) 



The maximum deviation from Model l is 71,704 pascals per meter, while the maximum deviadon 
from iModel II of the data shown in Figure 8 is 40,258 pascals per meter. In an effort to understand 
if the variation can be attributed to any of die recorded test parameters, an analysis of variance 
(ANOVA) is conducted for air velocit>', material diameter, material densit>% ppij and material 
thickness. The following hypothesis was tested: 



Ho: There is no difference in the model deviation means 

Ha: At least two of the model deviation means differ. 

The complete results from die ANOVA along with die histograms of the variable categories are in 
Appendix III. The analysis shows diat there is sufficient evidence (at 99% confidence) to reject Hb 
for all die variables categories investigated in Model 1. Rejection of the null hvpodiesis indicates that 
die model deviation average for die categories differ; dierefore some of die model deviation can be 
attributed to each of die categories investigated. The ANOVA for Model II indicates that the null 
hypothesis cannot be rejected for density and thickness, suggesting diere is not sufficient evidence 
that die model deviation can be attnbuted to diese two parameters. The analysis is summarized in 
Table 1 . 



Tablel: ANOVA Results 



Recorded Parameter 


Model I 


Model II 




F Statistic 


Conclusion 


F Statistic 


Conclusion 


Density 


4.785 


Reject Ho w/ 99% Conf. 


0.338 


Accept Ho 


PPI 


3.38 


Reject Ho 95% Conf. 


5.696 


Reject Ho w/ 99% Conf 


Air Velocity 


20.075 


Reject Ho w/99% Conf. 


6.70 


Reject Ho w/ 99% Conf 


Thickness 


6.3526 


Reject Ho \v/99% Conf. 


2.081 


Accept Ho 


Diameter 


12.382 


Reject Ho w/99% Conf. 


5.113 


Reject Ho w/ 99% Conf 



Based on the ANOVA results, an adaptadon of Model II is developed in an effort to capture as 
much of die variance in die data as possible. iModel III accounts for die variadon in hydraulic 
diameter, allows die velocity exponents to vary, and also allows for dissimilar exponents in PPI and 
relative density. The model coefficient of determination is approximately 99% and takes die form: 
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(18) 



ii = 3.176x10-^ p.^r^"'"^"'"'" 
K 



(19) 



Cp = 4J64xlQ-'p'J'ppi''\ 



(20) 



Eq. 18 indicates that as tlie diameter is decreased, the pressure drop increases. The diameter 
interaction seems to make physical sense as well as helping to increase the coefficient of 
determination for the model. The graph of Model III widi die experimental results is shown in 
Figure 10. 



Figure 10: Model III Pressure Drop vs Experimental 
Pressure Drop 
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Scrutiny of 'Eq. 9 

The next step in the analysis is to determine the feasibility of using Eq. 9 to determine die pressure 
drop dependence. Figure 11 shows die relationship between die pressure drop squared as a function 



9 



of velocity. The cun^e does not fit the physics of the problem and allows a negative pressure drop at 
low flows. The same trend is observed at foams wdth higher ppi number and reladve densit^^ these 
charts can be found in Appendix IV. For this reason, Eq. 9 is excluded from further scruuny and 
SHOULD NOT be used to calculate the pressure drop in metal foams. The chart also illustrates that 
as the pressure drop increases, Eq. 9 does a good job at determining the velocirvVpressure drop 
relationship. Eq. 9 was developed for use in sintered metal beds radier dian metal foams. Since the 
pressure drop in sintered beds is much higher than in metal foams die equadon probably works 
much better in that applicadon. 



Figure 11: Pressure Orop^ as a Function of Inlet Velocity 
5 ppi. Various Relative Densities 




Scrutiny of Eq. 1 1 

The last reladonship studied is Eq. 11. After looking at Eq. 18, one can see that the most dominant 
term in the pressure drop equation is U^. Therefore, it is not a far stretch of the imagination to 
believe that a single velocity term could yield a good fit to the experimental data. Figure 12 shows 
the experimental data fit with a power law curve. 



Figure 12: Pressue Drop as a Function of U*^ 
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Another interesting relationship can be derived from the data in conjunction with another 
experimental measurement Stainless steel metal foam was photographed in order to obtain the 
dependence of ppi and density on the strut diameter. A grid was over-laid on the sample, the pore 
size and strut diameter was measured with a ruler and scaled appropriately. Figure 13 shows a 
picture of the cell and the connecting struts. 



Figure 13: Microscope Picture of Metal Foam 
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Based on the results, the strut diameter is proportional to relative density^ and inversely propordonal 
to ppi. The relationship is found to follow, ' 



W - 1 A Pret 

^"^•^ (21) 

ppi ^ ^ 



Figure 13 shows a good fit widi the experimental data. 

Figurel3: Strut Diameter Model Compared to Actual Data 




The pore diameter is related to die strut diameter as: 




The comparison of die actual to calculated in shown in Figure 14. 



Figure 14: Experimental and Calculated Pore Size as aa Function of Strut Diameter 
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Finally, based on Eq. 10, Eq. 11, Eq. 21, and Eq 22, we derive the dependence of pressure drop on 
ppi, relative density, and velocity. The equation takes the form: 



0.4-0.1/n 
Pre I 



a r^a 



L(l-a) 



2-m 



rr2-m/ _,-0.7+0.7m'\ 



(23) 



Using a curve fitting algorithm, we see that m-0.19 and | -182.595. 



Comparison of Stain/ess Steel Pmsun Drop to FeCrAlY Pressure Drop 

Since the sintering properties of stainless steel and FeCrAlY are similar, one would expect that the 
strut diameter, pore size, and consequendy, the pressure drop, would be similar as well. The strut 
diameter and pore size was determined by Cambridge Univerist)' for FeCrAlY using image analysi 
The dependence of the strut diameter on ppi and reladve density from their report, is found to be 



^3 



= 0.601359-^. (24) 
ppi ' 



Eq. 24 is subtracted from Eq. 21 and dien divided by Eq. 21 to give the reladve difference in strut 
diameter and is shown in Figure 1 5. The plot indicates diat die relative difference is small at low ppi, 

but increases as ppi increases. 

Figure IS: Relative Difference in Strut Figure 16: Reladve Difference in CeU 

Diameter Compared to ppi Diameter Compared to ppi 

and Relative Density and Relative Density 




Also, the cell diameter is found to fit the equation, 



a = 



1.18 



Pre/ 

100 
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(25) 



Eq. 25 is subtracted from Eq. 22 and dien divided by Eq. 22 to give the reladve difference in cell 
diameter and is shown in Figure 16. Again, the plot indicates that the reladve difference is smaU at 
low ppi, but increases as ppi increases. 

Due to the calculated differences in the strut diameter and cell size, Eq.lO is used in conjunction with 
Eq. 24 and Eq. 25. The resultant equation, 



-jO.31-0.19m 
l_g0.04 



g r'a 



V)- 



0.4l+0.4lm 



(26) 



is used in a curve fitting algoridim indicating ^ = 4.219x10* and W = 0 . Aldiough the values of 
die constants in Eq. 26 and Eq. 23 are much different, the results from the equations are not much 
different for one anodier. Figure 17 shows die results of using die stainless steel correlation (Eq. 23) 
to predict the pressure drop in FeCrAlY. 



Figure 17: Using the Stainless Steel Calculation to Predict the Pressure Drop in 

FeCrAlY 
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The chart shows that the stainless steel prediction posses the same trend and magnitude as the 
FeCrAlY prediction. Since the pressure drop of die two materials is so much alike, one would expect 
that the material pore sizes and cell diameters are much alike as well. Therefore, the data suggests 
that the strut diameter and pore size measured manually differs from d:iose obtained from image 
analysis. Further measurements should be taken to understand why there is a difference. 

Experimental Results For Differing Inlet and Exit Diameter 

In all the preceding results, the inlet diameter and exit diameter of the part was kept constant diereby 
keeping the air velocity constant The question arises, "If the velocity of the stream through a metal 
foam changes along its length, which velocity should be used in the above equations, the inlet, the 
exit, or an average?" In order to answer this question, an apparatus was constructed widi an inlet 
orifice placed in front of the metal foam part, shown schematically in Figure 18. 

Figure 18: Schematic of Testing Apparatus 

4 4 4 4 ^ vout 

Vin 

The inlet orifice diameter, inlet velocity, metal foam diameter, metal foam diickness, metal foam 
relative density and ppi were all varied. The experimental pressure drop versus die calculated 
pressure drop was plotted using die inlet velocity and exit velocity. These plots are shown in Figures 
18 and 19. 



Figure 18: Experimental Pressure Drop 
Compared to Calculated w/ Inl t 
Velocity 
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One can see from Figure 18 that the inlet velocity gives fair accuracy at low velocities, but at high 
inlet velocities the pressure drop is overestimated. 



Figure 19: Experimental Pressure Drop 
Compared to Calculated w/ Exit Velocity 
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Conversely, one can see that the pressure drop is underestimated when using the exit velocit}- in the 
pressure drop equation. Due to the inaccuracy in the pressure drop prediction, it becomes necessary 
to develop a velocit)'-averaging scheme to generate a velocity to be used in the pressure drop 
equations. 



As an initial pass, a simple average is taken CUfin+Utbui)/2. The results are shown in Figure 20. 
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Figure 20: Experimental Pressure Drop Compared to 
Calculat dw/ Average V locity 



1600000 
1400000 
1200000 

? 

2 1000000 
a 

Q 800000 
I— 

8 600000 

Q. 

400000 
200000 
0 







♦ Experimental Data 


D 
♦ 




□ Calculated w/ Average 
Velocity 


g 

□ 


□ 








□ 


□ 




♦ 




a 


R 








"S — 

i 


9 













10 



20 30 40 

Inlet Velocity (m/s) 



50 



60 



70 



In order to determine if the error can be improved on, the weighted average function, 
along with the weighting function. 



k = 



1- 



200 



-1.8, 



1- 



Prel 

100 



1-- 

1 



0.74 I- 



(27) 



(28) 



is used to calculate the pressure drop. The values for the exponents were arrived at through a curve 
fitting algorithm which maximized the value when the calculated pressure drop was plotted with 
the experimental pressure drop. These results from the more complex averaging function are shown 
in Figure 21. The chart indicates that the simple averaging function is not a bad approximation when 
determining the pressure drop. 



Figur 21: Experim ntal Pressure Drop Compared to 
Calculated w/ Weighted Av rage V loclty 
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VELOCITY PROFILES IN IRON BASED METAL FOAMS 

In applications where metal foams are used for flow diffusion, die velocity profiles at the exit of the 
material are of interest. Measurements were taken in an attempt to determine the exit v'elocity profde 
and its dependence on the incoming velocit}^ the part diameter, part length, the ppi, and the relative 
densit)^ Figure 22 and 23 show a typical velocity profile for a 5ppi, 12% relative density' stainless 
steel part with a '/z inch inlet orifice and a 3.75 in part exit. If the velocity is plotted as a function of 
distance, a Gaussian profile emerges (Figure 24). The maximum velocit)'^ and the velocity' scatter can 
then quantify the velocity contour allowing us to see trends when parameters are varied. In general, 
the profile distribution is assumed to follow 

U, = Ur (29) 
The values of U j- and 2a are shown in Figure 24 for the particular case. 
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Figure 22: Exit Velocity Contour for a 5ppi, 12% Relative Density Stainless Steel w/ 0.5" 

Inlet Diameter and 3.75" Exit Diameter 




Figure 24: Exit Velocity Profile for a 5ppi, 12% Relative Density Stainless Steel w/ 0.5" Inlet 

Diameter and 3.75" Exit Diameter 
Average Velocity (Blue Stars) Compared to Gaussian Distribution (Green Line) 
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We can determine U j-^ by conservation of momentum and by integrating die proposed velocity 
profile. U is shown below in terms of the other variables. 



Ur =2 



D 



(30) 



It also reasonable to believe that 2a is a function of the foam thickness, ppi, relative densit)', inlet 
velocity and diameter. A multiple variable non-linear curve fit was conducted to determine the 
dependence of 2a on die independent variables. The following relationship emerges. 



2ct =3.513 



D, 



1 



rrOM __;n.6 -.0.66 I. 
^/,. PPl Prel ^ 



0.6 ^0.66 r3.09 



(31) 



however, die coefficient of determination is small (R2=30%) and can yield significant scatter from 
die experimental data. It is recommended diat computational fluid dynamics be used if accurac}' is 
required. Eq. 31 is intended to give a die flow scatter trends radier dian die acmal values. 
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SUMMARY 



The relationships for pressure drop as a function of velocity are presented. The set of equations with 
the largest coefficient of determination at 99% is: 



AP _ 1 



^U'+pCU' 
K 



(18) 



ii = 3.176xl0-^p^,r^"'"^'''"°' 
K 



(19) 
(20) 



Cp = 4.764xl0->;;>p/''\ 

Eq. (18), (19), and (20) yield good results over a range of PPI, relative density, thickness, and part 
diameter. In general Eq. (19) and (20) vary depending on die material. Therefore die values of — 
and Cp must be regressed for the different materials used. 

With some loss in accuracy, but gaining a generalized formula for all materials, die following equation 
is used: 



^p . 1 



a r^a 



,{l-a) 



l2-ot 



(10) 



where the cell diameter and the strut size of die material are the variables used to dictate die pressure 
drop. 

The strut size of stainless steel id found to follow the relationship, 



d = \A 



ppi 



0.7 



(21) 



and the cell diameter can be described by 



a = 




l_e0.04 



(22) 



There are conflicting results for die stmt si2e and cell diameter of FeCrAlY when calculated manually 
and using visual analysis software, but using Eq, (21) and (22) in Eq. (10) yields approximately die 
same results as the experimental data for die pressure drop. Therefore, it is assumed diat die strut 
diameter and cell diameter are close to those calculated for stainless steel. 

In die case where die inlet and exit velocities differ, a velocity averaging scheme should be used. The 
best fit to the experimental data was found widi the weighted average, 



(27) 



and the weighting function, 



k= 1- 



PPi 
200 



-'•V 



J Prel 
100 



V 



1-- 
1 



(28) 



\ ' J 

The velocity proBe at die exit of the metal foam is found to follow the Gaussian distribution, 



Jr JmMX 



(29) 



The constant Ufm„ is found through conservadon of mass and die integration of Eq. (29) as 



U. =21 

-/mu 



u 



(30) 



2a is a function of die inlet velocity, die inlet diameter, die exit diameter, die PPI number, die 
relative density, and die diickness. A multiple variable non-linear curve fit yields die following 
relationship. 



2(T=3.513f:^ 



vi.08 



1 



^^:0-6 ^0.66 r3.09 



(31) . 



The R2 value for this curve fit was low, but gives the tends for die maximum velocity and flow 
dispersion within metal foams. 

Furdier testing on odier materials such as titanium and nickel is being carried out to fully develop 
pressure drop and flow profile equations for a wide range of alloys. 
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Appendix I: Pressure Drop as a Function of Velocity 



Appendix II: Coefficients for Velocity Terms in Eq. 4 



rrl iNumOCT 


Relative Density (%) 


M/K (N s/m') 


pC (kg/m") 


5 


3.97 


160.7396 


l.OOE-OS 


5 


4.66 


221.7 


218JS 


5 


9.2 


349.71 


260 


10 


4.18 


194.74 


l.OOE-05 


10 


4.31 


242.63 


420 


10 


16.5 


654.86 


731.66 


20 


4.7 


286.23 


392.62 


20 


7.7 


739.84 


619.62 


20 


9.5 


442.16 


75.754 


20 


10 


607.84 


1451.5 


20 


12.5 


944.39 


694.72 


30 


4.1 


342.94 


96.617 


30 


4.4 


431.82 


180.24 


40 


4.8 


482.94 


0 


40 


5.5 


508.61 


715.18 


50 


4.5 


557.02 


870.11 


60 




585.67 


841.57 


60 


5.4 


720.62 


2257.4 


60 


7.9 


565.38 


5287.9 


60 


8.9 


1969.2 


2437.8 


60 


11.1 


2194.1 


410.47 


60 


13.3 


11220 


4699.2 


60 


16.7 


110.97 


7722 


60 


7.1 


1910.4 


3505.1 
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Appendix III: Avova Analysis 



Appendix IV: Pressure Drop Squared as a Function of Velocity 
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